Breaking of vortex lines - a new mechanism of collapse in hydrodynamics 
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A new mechanism of the coUapse in hydrodynamics is 
suggested, due to breaking of continuously distributed vor- 
tex lines. Collapse results in formation of the point singu- 
larities of the vorticity field |i7|. At the collapse point, the 
value of the vorticity blows up as {to — t)~^ where to is a 
collapse time. The spatial structure of the collapsing distri- 
bution approaches a pancake form: contraction occurs by the 
law h ~ (io — i)^ along the "soft" direction, the characteris- 
tic scales vanish like h ~ (io — t)^" along two other ("hard") 
directions. This scenario of the collapse is shown to take place 
in the integrable three-dimensional hydrodynamics with the 
Hamiltonian Ti, — J \fl\dr. Most numerical studies of col- 
lapse in the Euler equation are in a good agreement with the 
proposed theory. 

PACS: 47.15.Ki, 47.32.Cc 



I. INTRODUCTION 

Collapse in hydrodynamics of an ideal incompressible fluid, 
as a process of singularity formation in a finite time, is one 
of the central problems in the theory of developed hydrody- 
namic turbulence. The classical examples of such type spec- 
tra are the Phillips spectrum for water-wind waves H and the 
Kadomtsev-Petviashvili spectrum for acoustic turbulence H. 
In the first case white caps - wages of water surface - play a 
role of singularities, and in the second case these are density 
breaks (or shocks). 

The question about collapse in hydrodynamics is an old 
problem. For example, in 1981 P.Saffman M] considered col- 
lapse as one of the most important problems in hydrodynam- 
ics (see also papers |g] and references therein), probably al- 
ready L.Richardson and A.N.Kolmogorov understood an im- 
portance of this problem. In spite of so long a history of 
the question, there is no deep understanding of the nature 
of a collapse in hydrodynamics, even though there are many 
numerical simulations testifying to the collapse existence col- 
lected by now. As far as theory is concerned, essential results 
are absent, and, moreover, there is no common agreement 
about collapse as a subject for incompressible hydrodynamics 
(see, e.g.. Sec. 7.8 of the book of U.Frisch pi and references 
therein). In the theory, the only appreciable exclusion is the 
work of V.E.Zakharov, 1988 |7| (more detailed publication has 
appeared in 1999 M), where the consistent theory of collapse 
for two anti-parallel vortex filaments of small thickness was 
developed in quasi-two-dimensional approximation, when a 
flow is almost two-dimensional with a slow dependence with 
respect to the third coordinate (see also the paper M). A 
significant progress in studying hydrodynamic collapse has 
been achieved in numerical simulations of the Euler equation. 
Many numerical experiments testify that the value of the vor- 
ticity \Cl\ becomes infinite in isolated points in a finite time. 
As follows from the papers of Kerr ml, Grauer, Marliani 



and Germashewsky |lj, Polz ||12[, Boratav and Pelz |13|, |r2j 
grows at the collapse point like {to~t)~^ , where to is a collapse 
time. According to |lO|, [h3|, a spatial scale of the collaps- 
ing distribution contracts as (to — t)^' ^. In the recent paper 
by Kerr h^ , an anisotropy of the collapsing region has been 
reported. The data processing gave two scales, one of them 
being contracted as the root h ~ (to — t)^'^, and another one 
as /2 ~ to — t. It should be noted that the initial flow either 
possessed a definite symmetry or it was close to a symmetric 
flow in most numerical simulations. As a result, several sin- 
gularities arose simultaneously. For example, the evolution of 
two anti-parallel vortex tubes was studied in jlOl. The col- 
lapse here is caused by the Crow instability |l4|] leading at the 
nonlinear stage to the vortex reconnection. Symmetry of the 
flow makes the collapse to happen in two symmetric points. 

In the present paper, we suggest a new mechanism of the 
singularity formation connected with breaking of continuously 
distributed vortex lines. This mechanism is not related to any 
symmetry of the initial vorticity. The collapse itself is possible 
in one separate point. Probably, just this type of c ollap se has 
been observed in the recent numerical experiment |15| . 

The mechanism suggested can be naturally incorporated 
into the classical catastrophe theory llq. From this point 
of view, collapse can be considered as caustic formation for a 
solenoidal field. It is not so easy to understand, how a collapse 
arises in the Eulerian description. First, this is connected with 
a hidden symmetry of the Euler equation, i.e^he relabeling 
symmetry (for more details see the reviews |l8|, |19|). This 
symmetry generates the conservation law for the Lagrangian 
invariant - the so-called Cauchy invariant, which is expressed 
through the velocity curl and the Jacoby matrix of a map- 
ping from the Eulerian variables to the Lagrangian ones, and 
by this reason this invariant occurs very nonlocal in terms of 
the velocity field. On the one hand, the Cauchy invariant is 
known as invariant which characterizes the property of frozen- 
ness of vortex lines into a fiuid. On the other hand, all known 
conservation laws for vorticity, such as the Kelvin's and Er- 
tel's theorems, conservation of the topological Hopf invariant, 
being a measure of the fiow knottiness, are a simple conse- 
quence of the Cauchy invariant constancy. The frozenness of 
vortex lines means that fluid particles are pasted to a given 
vortex line and never leave it. A destruction of frozenness 
is possible only due to viscosity, i.e., beyond ideal hydrody- 
namics. Therefore, as a next natural step in the vortex mo- 
tion description, a mixed Lagrangian-Eulerian description has 
been introduced, where the main object is a vortex line [ pOJ , 
pi[ . Each vortex line in this description is labeled by a two- 
dimensional Lagrangian marker, while the third coordinate 
serves as a parameter determining the curve. This represen- 
tation, which we called as the vortex line representation, is a 
key-point in the description of hydrodynamic collapse which 
can be considered as a process of a caustic formation for the 
solenoidal vorticity fleld. 

The paper is organized as follows: in Sec. II we introduce 



the vortex line representation and explain its meaning. In 
Sec. Ill we consider the 3D integrable hydrodynamic model 
introduced in our previous paper [po|. The Hamiltonian of 
this model is unusual, it is expressed through the absolute 
value of the vorticity \ft\ 



Hamiltonian coincides with the kinetic energy of the fluid 



H 



\n\dr. 



fl.l) 



The given model can be integrated by means of combina- 
tion of the vortex line representation and the inverse scat- 
tering transform. By applying the vortex line representa- 
tion, the Hamiltonian is decomposed into a sum of Hamiltoni- 
ans of non-interacting vortex lines. Dynamics of each vortex 
line is described by the integrable one-dimensional Landau- 
Lifshitz equation for a Heisenberg ferromagnet or by its gauge- 
equivalent - the nonlinear Schroedinger equation. Thus, the 
integrable hydrodynamics represents a hydrodynamics of free 
vortex lines. As for hydrodynamics of free particles - hy- 
drodynamics of dust with a null pressure (see, e.g. [E2|), for 
the 3D integrable hydrodynamics typical singularities are also 
caustics. For hydrodynamics of dust, density turns into infin- 
ity at caustics. Unlike the dust density, being a scalar charac- 
teristics, vorticity is a vector solinoidal field. Therefore, the 
latter imposes some restrictions on a spatial structure near 
singularity. As it is shown in Sec. IV, the singularity struc- 
ture, being very anisotropic, turns into a pancake form. The 
spatial collapsing distribution at t -^ io leads to quasi-two- 
dimensional. Along the "soft" direction a more rapid com- 
pression takes place as h ~ (io — t)^ , and along two other- 
"hard" - directions h ~ {to — t)^". At the collapse time, 
the vorticity vector lies in the pancake plane and its value 
ri| blows up like (io — t)~^. This behavior corresponds to 
the general situation. The degenerated case is considered in 
Sect. V, where we consider collapse for topologically nontriv- 
ial axi-symmetric distribution of vorticity in the form of the 
so-called Hopf mapping when any two vortex lines are linked 
once with each other. In this case two eigen-values of the Ja- 
cob! matrix vanish simultaneously in the collapse point. As 
its sequence, vorticity occurs to have more strong singularity: 
|f2| ~ (to — t)~'^. And we conclude in Sec. VI with discussion 
of numerical experiments on collapse observation in the Euler 
equation and their correspondence to the proposed theory. 



II. VORTEX LINES REPRESENTATION OF 
HYDRODYNAMICS 

Let us consider the equations of the hydrodynamic type 



'dt 



curl 



curl—— X Si 



(2.1) 



where 7i{i7} is the Hamiltonian of a system, fl{r,t) — 
curlp(r, t) is the generalized vorticity, p is the canonical mo- 
mentum. The vector field 



V = curl{SH/Sn) 



(2.2) 



is nothing else but the fluid velocity. By the definition 
div V = 0, i.e., we deal with an incompressible fluid. If the 



n^ / Y* = 



n(ri)-n(r2) 
Syrlri — r2| 



dridr2 



then the expression (2.2) yields the usual relation fJ = curlv 



between velocity v and vorticity CI, while Eq.(2.1) transforms 
into the Euler equation for vorticity 



'dt 



curl[v X n], divv = 0. 



The important property of the equation ( |2.1[ ) is the frozenness 
of vorticity into the substance, i.e. all Lagrangian fluid par- 
ticles at i > remain at their own vortex line. After that, it 
is natural to introduce a mixed Lagrangian-Eulerian descrip- 
tion when each vortex line is labeled by its own Lagrangian 
marker u, which lies in a fixed two-dimensional manifold A/", 
while a parameter s along the line has a meaning of an Eu- 
lerian variable. In such vortex line representation vorticity is 
expressed as follows |20| 



n{r,t) = / D,o{iy)du 
JaT 



5{r-Il{s,iy,t))^ds. (2.3) 



Here, the closed curve r — R(s, i^,t) corresponds to each vor- 
tex line v, so that Rs is its tangent vector. The quantity 
iloiv), with Rs being its tangent vector. The fixed function 
r2o(i^) is the strength of vortex loop. However, without loss of 
generality, this function can be put equal to the unity. This 
can be achieved by both an appropriate re-definition of labels 
V and changing the vortex orientation to the opposite one for 
those lines from the manifold A/", for which f2o(j^) < 0. There- 
fore, in the next section we will omit the multiplier Qo{i^) in 
front of d^u in the correspondi ng fo rmulae. 

The generalization of Eq. (2.3) to the case of arbitrary 
topology of the vortex lines is served by the formula: 

n(r, t)= I 5(r - R(a, t))(no(a) Va)R(a, t)d'a, (2.4) 

where J7o(a) is the Cauchy invariant, characterizing the 
frozenness property. The divar2o(a) ~ condition guar- 
antees automatically incompressibility for the field fl{r,t): 
divf2(r) = 0. In the expression (2.4), the vector 



b = (no(a)V„)R(a,i) 
is a tangent vector to the vortex line at the point 

r = R(a,t). (2.5) 



In the representation (2.3), an arc-length of line of the initial 
field i7o(a) can serve as the parameter s. 

After integrating (2.4) over a-variables, the vector f2(r, t) 
is expressed through the Jacobian J of the mapping (2.5) 
J — det||9R/9ajj and the Cauchy invariant f2o(a): 



n(R) = -(no(a)Va)R(a) 



(2.6) 



It is important to emphasize, that in this expression the Ja- 
cobian is not to be equal to the unity: J 7^ 1. Nevertheless, it 
does not contradict to the incompressibility condition for the 
fluid. 



As it was shown in [EOlpiJ, the equations of motion for 
vortex lines can be obtained directly from the equation of 
frozenness (^.1) 



[{(no(a)VJR(a,i)} x {Rt(a, t) - v(R(a, t), t)}] = 0. (2.7) 

This equation describes the transverse dynamics of vortex 
line: obviously any motion along a curve leaves the curve un- 
changed. In particular, this helps to understand why there are 
no restrictions imposed on the value of the Jacobian J. Let 
us recall that in a purely Lagrangian description of an incom- 
pressible fluid Jacobian is eq ual t o the unity identically. But 
accordingly to the equation (2.7), the motion of Lagrangian 



particles along vortex lines is excluded from the mapping ( [2.5| ) 
for the mixed Lagrangian-Eulerian description. It is exactly 
the reason why the Jacobian is not necessarily equal to unity 
now. This point is principal and will be instrumental below 
in explai ning how collapse is possible in the hydrodynamic 
systems (2.1). 

As it was shown in pd,El|, Eq.(2.7) can be written in the 
Hamiltonian form 



[{(no(a)Va)R(a,t)}xRt(a) 



SH{n{Il}} 



(5R(a) 



(2.8) 



This equation describes a motion of vortex lines in systems 
with an arbitrary Hamiltonian that depends on R through 
the i7(r, t) only. 

It is useful also to keep in mind that the expressions for 
such important characteristics of the system as its momentum 
P = fpdr and angular momentum M = J[r x p]dr, being 
transformed by integration by parts to a form, where, instead 
of p, the vorticity J7 is employed, and being then rewritten in 
terms of vortex lines, have the form 



P ~ - / [r X fl]dr = 



M ~ - / [r X [r X fl]]dr - 



^2 1 

a V — 



M 



[R X Rs]ds 



(2.9) 



u 



d^v - I [R X [R X RJ]ds 



(2.10) 



The ~-sign in these relations means that equalities take place 
up to integrals over surface with infinitely large radius. Hence 
one can see that the momentum and the angular momentum 
are composed of momenta and angular momenta of each vor- 
tex line, the momentum of a closed line being equal to the 
oriented area of a surface tightened on the vortex loop. 

It is easily to verify that uniform shift Ro of R does not 
change the momentum, while the angular momentum is sub- 
jected to the well known transformation 



M 



M' = [Ro X PI + M. 



(2.11) 



model was introduced in our previous paper [E0| . The Hamil- 
tonian of this model is expressed through the absolute value 
of n(r,i) 



H= / |n(r)|dr. 



(3.1) 



and the equa tion of motion coincides with the frozenness 
equation (2.1) with velocity 



V = curlr 

where r = (ft/Q.) is the unit tangent vector along the vortex 
line. Assuming all the lines closed, choosing the labeling by 
such a way so tha t Q o(u) — 1, and substituting the represen- 
tation (2.3) into (3.1), it is easy to see that the Hamiltonian 
is decomposed as a sum of Hamiltonians for the vortex lines 



H{R.} 



dR 



ds 



ds. 



(3.2) 



Here, the integral over s is the length of the vortex line with 
the index v. The equation of m otion for the vector R(!/, s), 
in accordance with the Eq. (2.8), is local in these variables - 



it doesn't contain an interaction with other vortices: 



[R. X Rt] 



T X T X Ta 



(3.3) 



By this reason, not only the total energy, momentum, and 
angular momentum are conserved, but also the corresponding 
geometrical invariants for each vortex loop: its length 



H{u) 



|Rs(!^)|ds, 



the oriented area spanned on the vortex loop which coincides 
with its momentumits momentum 



[R(!/) X 'R.s{u)]ds, 



and its angular momentum 



1 



M(i/) = - [R(i/) X [R(i.) X Yis{i^)]]ds. 

It is important to pay attention to the following fact: The 
equation ( p.3| ) is invariant with respect to changes s — > s{s, t). 
Therefore, it can be solved for Rt up to a shift along the vortex 
line - the transformation leaves the vorticity CI unchanged. 
This means that to find the vorticity fJ it is enough to have 
one solution of the equation 



IRslRt = [TXTs]+l3Yis 



(3.4) 



which follows from Eq.(3.3) for some choice /3. This leads 
to an equation for t as a function of the filament length I 



III. INTEGRABLE HYDRODYNAMICS 

In this and two next sections, we will show how and why 
collapse is possible in 3D integrable hydrodynamics. This 



^It is worth to notice that this property is common 
for all systems with the Hamiltonians of the type Ti — 
J F{t, {t\7)t, {rVfr, . . .)|n|dr. To explain the idea of col- 
lapse of vortex lines, we have chosen the simplest example 
(|3.l|), which has a physical meaning. 



{dl = |Rs|ds) and time t (by choosing a new value f3 — 0), 
which reduces to the integrable one-dimensional (ID) Landau- 
Lifshits equation for a Heisenberg ferromagnet 



dr 
dt 



dP 



(3.5) 



This equation, in its turn, is gauge equivalent to the ID non- 
linear Shrodinger equation |23[ 



ii^t+i'ii + {l/2)\i,fij = 



(3.6) 



and, for instance, can be reduced to the NLSE by means of 
the Hasimoto transformation fE4| 



i / xil,t)dl 



ip{l,t) = K{l,t) ■ exp 



where k(Z, t) is a curvature and x{ly *) the line torsion. 

The system under consideration has direct relation to hy- 
drodynamics. As it is known |25|], pq], the local induction 
approximation for a thin vortex filament, under assumption 
of smallness of the filament width to the characteristic longi- 
tudinal scale, leads to the Hamiltonian (3.2), but only for a 
single separate line. The essence of this approximation is in in 
replacing the logarithmic interaction law by a delta-functional 
one. When the widths of the filaments are small comparable 
with distances between them, in the same approximation, the 
Hamiltonian of vortex lines transforms into the sum of the 
Hamiltonians of independent vortex loops, yielding in a "con- 
tinuous" limit the Hamiltonian (5.1). 

By such a way, we have the model of 3D integrable hydro- 
dynamics of free vortex filaments. In this model, each vortex 
is a nonlinear object with its own internal dynamics. As we 
will see later, already in the framework of this model, a singu- 
larity formation is possible. Singularities appear in this model 
as a result of intersection of vortex lines that is analogous to 
the phenomenon of wave breaking in gas-dynamics. 



A. Stationary vortices 



Let us consider now the simplest solution of Eq.(3.3), i.e., 
a stationary propagation of a closed vortex line: Rt = V = 
const. In this case the velocity V is determined from solution 
of the equation 



[Rs X V] = [f X [f X Ts 



(3.7) 



It is easily to check that this equation follows from the vari- 
ational principle 



5{H{u) -V-Piu))^0, 



(3.8) 



i.e., any solution of (3.7) represents a stationary point of the 
Hamiltonian for a fixed momentum P(!^). The equation (3.7) 
can be simply integrated, being rewritten in terms of the bi- 
normal b and the curvature k of the line as follows 



that gives 



[r X V] = k[t X b] 



Kb. 



(3.9) 



(3.10) 



A constant value of the velocity V in this expression implies 
constancy of the curvature k, i.e. the vortex line must be a 
ring of radius r — 1/k and 



V = l/r 



(3.11) 



The direction of the ring motion is perpendicular to its plane. 
It is interesting to note that the exact answer to the velocity of 
a thin (with width d <^ r) vorte x ring in ideal hydrodynamics 
( p^) coincides with Eq.( |3.1C| ) up to the logarithmic accu- 
racy that just differs the considering model from the Euler 
equation. 



Stationary solutions (3.10) in the form of rings are remark- 



able within this model, because they are stable, moreover, 
they are stable in the Lyapunov' sense. Remind, that momen- 
tum P of a closed vortex line is its oriented surface spanned 
on the loop: 

P^Sn, 

where S is the surface value, n its normal. Inasmuch as the 
Hamiltonian of a vortex loop coincides with its length, a max- 
imum of the momentum, or, that is the same as a maximum 
of surface S is obviously attained, for fixed length, at the per- 
f ect c ircle. Just this proves stability of the vortex ring solution 
.101) in the Lyapunov sense. 



IV. COLLAPSE 



The solution (3.10), (3.11) enables us to construct the sim- 
plest mappings R = R(z^, s,t). 

Let all vortex lines be circle-shaped and oriented in the 
same direction, for instance, along z-axis. We will see fur- 
ther that collapse in our model is a purely local phenomenon. 
Therefore, it is sufficient to consider some vortex tube (which 
can be imagined as a torus) to find a mapping. Let vortex 
rings be distributed continuously inside the tube. We label 
each vortex line by the two-dimensional parameter v, which 
values coincide with coordinate of some cross section of the 
tube at t = 0. We will use the ring arc-length as longitudinal 
parameter s {ds = rd(j>, where 4 ) is t he polar angle around 
z-axis). Then, with the help of ( 3.10 ), the desired mapping 
can be written as follows 

R — Ro(!^) + r{iy)cos<l)ex + r{h')sm(l)ey + V{iy)tez. (4.1) 

In this formula ex,y,z are unit vectors along the corresponding 
axes. 

It can be easily verified for this mapping that the Jacobian 
is a linear function of time 

J = y ' ^' ^,1 = M'^, s) + A{u, s)t. (4.2) 

Here A{i', s) is a coefficient linearly dependent on the veloc- 
ity derivatives with respect to i/ and Jo the initial value of 
Jacobian. 



Dependence J (4.2) on time means that for every fixed pair 



v and s there exists such a moment of time t — tiy, s) ( i > 0, 
or t < 0), when Jacobian is equal to zero: J{u,s,t) — 0. 
Denote as to the minimal value of i = i(i/, s) at i > 0. And 
let this minimum be attained at some point a = ao (here we 



denote a point {ui, 1^2,3) as a). It is evident that at t = to 
{dt/da) \a=ao = or 



9J(a, t) dt 



(4.3) 



V.J(a,t)|„^„^ + 



dt da 



It is clear also that at f = to the tensor of second derivatives 
of J against a, 

oaiOaj 
will be positive definite at the point a = ao. Hence, it is easy 
to define a behavior of the Jacobian in a small vicinity of 
a = ao . Expansion of J near this point (in a typical situation) 
at i ^ to is as follows 



J(a, t) = a{to —t)+ jijAaiAuj + . 



(4.4) 



where 



dJ{a,t) 



dt 



> 0, Aa = a — ao. 



t^tf) ,a^aQ 



Those are the leading contributions to the Jacobian expansion 

Geometrically, the above expansion corresponds to a suffi- 
ciently simple picture. The (hyper-) surface J = J(a, t) de- 
forms with time in such a manner, that its minimum reaches 
the (hyper-) plane J = at t = io, where two surfaces touch 
each other. Obviously, for smooth mappings in a typical case, 
this touching takes place in one separate point. In a degener- 
ated situation, touching is possible in a few points simultane- 
ously, or even at a curve. A case, when two eigenvalues of the 
Jacoby matrix J tend to zero simultaneously at the collapse 
point, will be regarded also as degenerated (such an exam- 
ple will be considered in the next section). In that case one 
should keep the next terms in the Jacobian expansion along 
the corresponding directions. We would like to repeat that 
all these cases can not be co nsid ered as typical ones. 

In accordance with Eq. (2.6), the equality J = at the 
singular point means the formation of a singularity for the 
vorticity at the moment t = to: 



n(r,t) = 



fio(t^)Rs 



a{to —t)+ -jijAaiAaj 



(4.5) 



It is important, that the numerator in this fraction - the 
tangent vector of a vortex line - does not vanish at the point 
ao due to its geometrical meaning. Therefore, the vorticity at 
the singular point blows up as (to— 1)~^, and the characteristic 
size of the collapsing distribution in a-coordinates decreases 
as V^o — t. 

The above type of collapse arises as a result of vortex line 
breaking when one vortex overtakes another. For flows of the 



^ For instance, the term containing mixed time-space deriva- 
tives, in the general case, is a small correction to the first term 



in (4.4) 



general type (without symmetries) a singularity must arise at 
the first time always in one separate point. 

As we will see further, for the given type of collapse, the 
dependence ( [4.q ) for fl{r,t), derived for the particular initial 
distribution, is actually the general answer, which can be ap- 
plied not only to the integrable hydrodynamics but also to 
the whole family (2.1) of hydrodynamic systems, of course, 
under the condition that they admit such (quasi-) inertial 
regime of collapse. What necessary conditions should satisfy 
the Hamiltonian of some particular system, in order to ex- 
hibit such regime? Now the answer to this question has been 
unknown that provides a wide field for future investigations, 
both theoretical and numerical. 



A. Non-stationary vortices 

Let us consider now the integrable hydrodynamics in a 
more general case when closed vortex lines are not circles. 
In this situation, it is possible to map each vortex contour to 
some vortex ring. It is then natural to introduce the (mean) 
direction n, as well as the mean area S = vrro, with the help 
of the expression for momentum of vortex line, P = nS. The 
position Ro of the ring center changes linearly in time. A cor- 
responding mean velocity of the motion of closed line must be 
directed along the momentum, in order to satisfy the conser- 
vation law for angular momentum. The mean velocity value 
Vo, generally speaking, is a function of those fundamental in- 
tegrals of motion, which are independent on the origin choice 
- the Hamiltonian (i.e. the length L), the momentum, and 
the projection on the vector n of the angular momentum. It is 
clear also that increasing the contour sizes in A times should 
result in decreasing the velocity in A times, so that 



V - s P 



u 



fWiT^P^ (M-P) 



L5 



(4.6) 



Explicit dependence of the U{^, ry)-function has been unknown 
now. Its first argument can be considered as a measure of a 
" crump leness" of the vortex line and changes in the limits 
< ^ < 1, and the second argument determines a measure of 
a spirality of the line. In a more or less reasonable approxima- 
tion, which doesn't lead to excessive errors, one can suppose 
that U{^,rj)^U{l,0) = 1. 

After introducing the mean characteristics, the mapping 
R(!/, s, t) can be represented as the sum: 

R(i^, s, t) = R(!/, s, t) + Sv{v, s, t), (4.7) 

where the mean value lA.{v, s, t) is given by the relation 



R(i/,s,t) = Ro 



- rocos (f)' e'x + rosin 0' e'y 



Ro = Vbn. 
(4.8) 



Here angular parameter 0' = 2tt{s/L) is proportional to the 
arc-length s, and unit vectors e'x,e'y lie in the plane, per- 
pendicular to the local z'-axis, directed along the n. The 
relationship between ro and Vo is given by the Eq. (4.6). 



The vector function Sr{v,s, t) describes deviations (generally 
speaking, not small) from the mean value R(!/, s,t). 

The separation of the mean and oscillatory motions, intro- 
duced by means of (p~7|), ([4.8[), (p~q) for each vortex contour. 



shows that the mapping H — I{.{v,s,t) at each fixed value 
a = {v, s) is a hnear function of time with nonlinear oscil- 
lations, which are described b y th e Landau-Lifshits equation 
( p.5| ) or its gauge equivalent (3.6). The linear mean depen- 
dence reflects the fact that the model under consideration is 
a model of free vortices. The collapse thus arises as a result 
of a running of one vortex into another. Due to a continuous 
distribution of vortex lines, their "density" increases infinitely 
in some point. 

A similar situation takes place in the model for large scale 
structure formation in the Universe studied by V.I. Arnold, 
Ya.B.Zeldovich and S.F.Shandarin [p2|. In the basement of 
this model, the suggestion lies about initial dust- like distribu- 
tion of masses, when their behavior can be described by the 
zero-pressure hydrodynamic equations 



pt +divpv = 0, 
0. 



- = V. + (vV)v 



(4.9) 
(4.10) 



The integration of this system in the Lagrangian variables 
gives that i) all fluid particles being free move with a constant 
velocity 



r = a + v(a)t. 



(4.11) 



and ii) the density p is express ed th rough the initial value and 
the Jacobian of the mapping (4.11) as follows 



P(r,i) = 



Po(a) 
J ■ 



(4.12) 



In the framework of this model, appearance of large scale 
structures is connected with the breaking phenomenon re- 
sulting in density singu larities, due to the Jacobian vanishing 
for the mapping (4.11). In a typical situation, these struc- 
tures have a pancake shape and can be considere d as pcurlo- 
galaxies. The formula, analogous to the Eq. (1.12), takes 
place, as we have seen above, for the vorticity $7 - s ee th e 
Eq. (2.6). However, there is a difference from Eq. (4.12), 
connected with the vector nature of the fJ fleld. 

In the given task we will be interested in a geometric struc- 
ture of a the singularity at i — *■ io, but t < io, i.e., in some 
sense, at the initial stage of collapse, but not at the developed 
stage, which has the meaning for astrophysics applications, 
but remains rather unclear for the incompressible hydrody- 
namics, when viscosity is necessary to be taken into account 
at small scales. 



is assumed to be non-degenerate: all its eigenvalues are posi- 
tive, and the matrix itself can be diagonalized. Hence it fol- 
lows immediately that compression along all principal axises 
in a-space will have the same law: la ~ \/io — i. Therefore, 
near the singular point, vorticity fi will have the self-similar 
asymptotics 



fi(r,i) 



no(!^)Rs 



(to -i)(a + 7«j'7«^j)' 



(4.13) 



where r\ — Aa/y/to — t are s elf-similar variables in a-space. 
However, the equation (4.13) does not mean, that compres- 



sion in r-space will be the same. 

When the Jacobian takes zero value, one of the eigenvalues 
of the Jacobi matrix becomes equal to zero. This eigenvalue 
(deno te it as Ai ), as it easily to see, coincides with the Jaco- 
bian (4^) at a small vicinity of the collapse point up to the 
almost constant multiplier A2 • A3. 

Represent now the Jacoby matrix J as a decomposition over 

An|^'"' >) and conjugated 
ems 



eigenvectors of direct {J\ip'-"^ >— ^ '"'''"' 



(< V;("'|J = A„ < ■(/;(") I) spectral proble 
_ dxk sr^ . ,(„) 7(„) 



(4.14) 



Here two sets of eigenvectors of direct and conjugated prob- 
lems are mutually orthogonal 

In a small vicinity of the collapse point eigenvalues A2,3 can 
be considered as constants, while 



Ai 



J 



A2A3 



= (A2A3) [a{to — t)+ "fijaiaj 



Here, for simplicity, we have placed the origin at the point 
a = ao. As for the eigenvectors, they also can be considered 
constant. 



Let us decompose the vectors x and Va in Eq. (4.14) 
through the corresponding bases, denoting their appropriate 
projections as X„ and A^. 



X„ =< xiV^ 



(n) 



>, 



d 
dAn 



<l/5'"'|Va> 



In this case the vector a is written in terms of An as follows 



E^ 



B. Structure of the singularity 

Let us consider in more details the structure of the col- 
lapsing region in a typical situation. First, it is clear from 
the above discussion that the vorticity distribu tion near sin- 
gularity will be given by the former expression (4.5). Second, 
the main features of the singularity will be determined by 
the Jacobian, that is the denominator of (4.5). Numerator 
(f2o(a)Va)R ( tangent vector to the lines) can be taken at 
the point a = ao, t = t p an d considered as a constant. 

According to Eq. (4.5) the Jacobian expansion contains 



positive deflnite symmetric matrix 7;^ of its second deriva- 
tives taken at the point a = ao, t = to- At t < to, this matrix 



As a result, Eq. ( 4.14 ) can be rewritten in the form 



3x2 _ 

dA2 ^ 



A2 



9^3 



A3 



(4.15) 
(4.16) 



Here, the matrix 

r„„ = 7«,4A2A3)-^^i"V^'"'l^'"'l'l^''"f , 

and parameter r — a(to — t)(A2A3)~^ is assumed to be small. 
It follows from here immediately that size reduction along the 
second {X2) and the third (X3) directions is the same as in 



the auxiliary a-space, i.e., ~ ^^ but along the "soft" direc- 
tion X\ behaves like r'^'^. Respectively, in terms of new self- 
similar variables ^i = X^Jt'-^I'', ^2 = X-^.jr^''^, ^3 = Xs/t"-^^, 
integration of the system gives for ^2 and ^3 a linear depen- 
dence on rj , while for ^1 it a cubic dependence 

^i = (l + r.,77,77,)r?i + irHry.77? + irn77?, i,i = 2,3 (4.17) 

C2 = A2r;2, Ca = ^sVs- (4.18) 



Together with ( [4.1^ ), the relations ( [4.17| ) and ( [4.18D deter- 
mine implicit dependence i7(r, i). The presence of two differ- 
ent self-similarities shows, that the spatial vorticity distribu- 
tion becomes strongly flattened in the first direction, taking 
a pancake form as t — > to- 

The direction of the field fl can be found from the incom- 
pressibility condition div $7 = 0. It is easy to see that in the 
leading order (as f — > io) the gradient of J is determined by 
the soft direction ei: 



VJ; 



-3/2 



dJ 



-ei. 



Contributions from two other directions are small in the pa- 
rameter r. 

Hence, it follows that vector lines of the field ft lie in the 
pancake plane that is in agreement with transversality of mo- 
tion of vortex lines (compare with Eg. ([2. 7])). 



V. EXAMPLE OF COLLAPSE IN THE 
DEGENERATED CASE 

In the previous section we have considered collapse for a 
non-degenerated situation, when the only one eigenvalue of 
the Jacoby matrix tends to zero at the touching point. Now 
we shall consider an example of collapse, when two eigenval- 
ues of the Jacoby matrix vanish simultaneously at the col- 
lapse point. We will examine an initial vorticity distribution 
with the nontrivial topology of vortex lines with linking num- 
ber A'^ = 1. This special distribution is the so-called Hopf 
mapping. There are several ways how to construct the cor- 
responding field ft. We will keep here the approach of the 
paper |2|. 

Following to that work, let us represent the field fl through 
the n-field (n^ = 1) 



n^{r) 



1 
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eQfj^(n- [dpn x d^n]), 



(5.1) 



where the n-field is supposed to be a smooth function of coor- 
dinates tending to the constant value e at the infinit y r — » 00. 
It is easy to check that in accordance with Eq.(5.1) each 
point n = no on the unit sphere S^ defines a closed vortex 
line . Indeed, parameterization of the unit vector n through 
the spherical angles 9 and ip allows to write down the field $7 
as follows 



n = —[Vip X VcosS]), 



(5.2) 



so that the variables if and cos 9 play the role of the Clebsch 
variables. Thus, each vortex line coincides in this case with 
the intersection of two surfaces ip = const and cos 9 = const, 



i.e. a closed vortex line is the pcurlotype in 7?.'' of a point 
on the sphere S^. By the definition, the Hopf degree N of 
mapping 7?.^ -^ S^ is called an integer number of linkages 
between arbitrary two vortex lines - pcurlotypes of two points 
on the unit sphere. The Hopf mapping (with N = 1) is given 
by the next relation: 



(n-a) = U{e-a)U\ U 



l-|-i(a- a) 



(5.3) 



where a are the Pauli matrices. 

Expressing the vector n from here and substituting the 
result into Eq.(p.l|), one can get (see p9|) 



r2o(a) = 



e(l 



2a(e ■ a) -|- 2[e x a] 



(l-ha2)3 



(5.4) 



As it was shown in Efl], all the fiux lines of this field are 
circles, and each line is linked once with another line. By this 
reason, the singularity formation is inevitable. 

It is worth to note that the field (5.4) has no singular points 
in the whole space and its absolute value depends only on the 
absolute value |a|. The unit vector r is defined everywhere: 



f(a) 



3(1- 



2a(e ■ a) + 2[e x a] 



(l + a2) 



(5.5) 



The velocity of each ring is connected with the binormal b(!/) 
and the radius r{iy) by the relation 

V(i/) ^h{u)/r{u). 

The radii r{v) of rings and their orientations b(i/) are integrals 
of motion. Only positions of centers of rings can change, and 
this motion occurs with a constant velocity for each ring. 

In the given problem, instead of the variables v and s, it 
is convenient to use directly the variables a, in which the 
mapping R(a, t) is written as 



R(a,i) = a-f tV(a) 



(5.6) 



where the velocity V(a) is expressed through the unit tangent 
vector r by means of the relation 

V(a) = [f(a) X (f(a)Va)r(a)]. (5.7) 

Hence one can find the expression for the Jacobian 



J(a, i) — det 



f fdV(a 



(5.8) 



It should be noted that the velocity of each ring is constant 
along the vortex line. Therefore the determinant of the matrix 
9V(a)/9a, which is the coefficient in front of f"^ in the previous 
expression, is equal to zero identically. As a result, J has a 
quadratic dependence on time t only: 

J = H-tei(a) +t^C2(a). 

Singularity formation corresponds to the case J ^ at some 
point. The collapse moment of time to will be determined by 
the coefficients ci and C2 so that to will be a minimal positive 
root of the equation 



minaJ(a, to) = Jmm(to) 



0. 



Calculations by means of Eq-s (5.5), (5.7) and (5.8) lead to 
the following expressions: 

V(a) = n^l-zyz ([^ X a](l - a^) + 2a(e ■ a) - 2ea^) , (5.9) 



J(a,t) = 



(1 + 0^)3 _ 8ia3(l + a^) + 4*2(1 + al - al - al) 



(l + a2)3 



(5.10) 



where as is the projection of the vector a to the axis e. 

Analyzing the latter expression one can show that at t < 1 
a minimum of the Jacobian is attained on the symmetry axis, 
i.e., at ai = 0, 02 = 0. In this case 



•J \axis 



(1 



+ 4(i-a3)2 



(l + ai)2 



(5.11) 



hence it becomes clear that the singularity takes place at to ~ 
1, as = 1, and the Jacobian tends to zero with the quadratic 
asymptotics. Thus, in this example 



inh 



{to-ty 



It can be easily verified also that the Jacoby matrix has two 
zero eigenvalues at the touching point which eigenvectors lie 
in the plane orthogonal to e. In the collapse vicinity the 
field n is directed along the vector e. Compression along this 
direction is linear in time I3 ~ (to—t), but in the perpendicular 
plane it is more fast, with the law h,2 ~ (to— i)^ • As a result, 
the singularity structure occurs strongly stretched along the 
anisotropy axis. 



CONCLUDING REMARKS 

Considering the hydrodynamic model with the Hamilto- 
nian Ti. = J \fl\dr, we have arrived at the conclusion that 
each vortex line in the given system moves independently of 
other lines. Just this property makes possible to form singu- 
larity in a finite time for the generalized vorticity fl(r, t) from 
smooth initial data. A typical singularity of this kind looks 
like an infinite condensation of the vortex lines near some 
point. Thus, the collapse in the integrable hydrodynamics 
has purely inertial origin. If one will assume that this type of 
collapse is possible also in the Euler hydrodynamics, then the 
asymptotics of the vorticity near the singularity point (the 
point of vortex "overturning") will be given by Eq-s (4.5) or 
(4.13) in a non-degenerated situation. That is already clear 



from the general consideration. Namely, the curl of the veloc- 
ity will blow up as {to — t)~^ . Exactly this dependence for 
vorticity near a singular point has been observed in practically 
all numerical simulations of the Euler equation, including the 
above cited: §, 0,0, 0, |i|, [@. 

However, not all the simulations are regarded to numerical 
integration of the Euler equations for continuous distributions 
of vorticity. The first numerical experiments |6| relate to in- 
vestigation of collapse for two anti-parallel vortex filaments 
which, as it was shown by Crow |l4], are linearly unstable 
with respect to transverse perturbations ( about the develop- 
ment of this direction see also |12]). The theory of collapse for 



two thin vortex filaments, as the nonlinear stage of the Crow 
instability, was developed by V.E.Zakharov [W], [p| (see also 
[b|). The conclusions of this theory are in a good agreement 
with the numerical experiments up to the distances compara- 
ble with a core size of filaments. This theory predicts decreas- 
ing distance between vortex filaments as -^to — t- For smaller 
distances the cores of vortices loose their round shape. They 
become fiat, and the process of attraction between vortices 
becomes more slow po|, |pl|. The same tendency was ob- 
served also in the numerical experiments of Kerr ||lO|, the 
most advanced (in our opinion) for the problem of recon- 
nection, and where, unlike |m], collapse of two anti-parallel 
vortices was simulated with continuously distributed vortic- 
ity. Besides a natural contraction of the minimal distance be- 
tween the distributed vortices, these simulations have shown, 
for the first time, the formation of two singularities in two 
symmetric points. While approaching a collapse moment in 
time, the explosive growth of the maximal value of vorticity 
was observed with the law {to — t)~^. According to the recent 
publications of Kerr ||l7| , the analysis of numerical data gave 
two distinguished scales, one scale being contracted as the 
square root: Zi ^ (io — i)^ , and another as the first power 
of time: I2 ^ to — t. 

In the work of Grauer, Marliani and Germaschewsky [|ll|, 
the successfuU attempt was undertaken to observe collapse 
for the initial condition not possessing a low symmetry. The 
initial vorticity was concentrated in the vicinity of a cylinder, 
and was modulated over the angle in such a way so that the 
simplest symmetries were absent. In the present time this 
experiment has the best spatio-temporal resolution. In this 
simulation appearance of a separate collapsing region was ob- 
served with the vorticity growth at the center as (to — t)~ ■ 

Thus, the results of the numerical simulations for ideal hy- 
drodynamics go along with the concept of collapse as a process 
of caustic formation for the solenoidal field. There is an agree- 
ment for behavior of the vorticity maximum. Concerning the 
spatial structure of the collapsing domain, we can make only 
a qualitative agreement. The results of the papers [M, |3ll] 
(where some contraction of the vortex core was observed at 
the initial stage of reconnection neglecting viscosity) and also 
the results of Kerr seem to support our theory. We would 
like to pay attention to the numerical results of the paper |33] 
where for short-time dynamics for the initial conditions in the 
form of the Taylor-Green vortex and for random initial con- 
ditions it was observed formation of thin vortex layers with 
high vorticity that also support our predictions. 

All these results allow one to say that the scenario pre- 
sented here looks like very plausible. 

We would like to repeat once more, that if such scenario 
takes place then behavi or o f th e vor ticity near of singular 
point is defined by Eq-s (4.5) or (4.13). The structure of this 
domain should be highly anisotropic: in one of two directions 
perpendicular to the vorticity, there is more rapid contraction 
(~ r'^") than with respect to other directions (~ r^"). The 
spatial distribution becomes close to the two-dimensional one 
familiar to a tangential discontinuity. The velocity of fiow in 
this region can be approximated with a good accuracy by the 
linear dependence 

V± ~ ^maxXl, 

i.e. the fiow looks like a shear fiow. The given type of col- 



lapse, according to the classification of |p2| , belongs to a weak 
coUapse: the energy captured into singularity (with account 
of the viscosity u) tends to zero as i^ — » 0. It is interesting to 
note that the dissipation rate ~ J Q'^dr from the collapsing 
region also vanishes when u ^ 0. 

One should note once more that unlike the model of free 
vortices considered here, in the Euler hydrodynamics vortex 
lines interact pairly in accordance with the Hamiltonian 



Ti-Euler — TT- 



lR(r.,s)-R(M,e)| '^ ^ 



(5.12) 



From the viewpoint of investigation of the collapse problem it 
is very important that the interaction function (i.e. the Green 
function for the Laplace operator) has the singularity for equal 
arguments R(/^, s) -^ R(Mj 0- If ^-his singularity would be ab- 
sent, i.e. if the interaction function would be completely reg- 
ular, then any initial vortex line distribution equivalent in the 
sense of iso-vorticity to a smooth field flo (a) , even very singu- 
lar distribution, would produce a sufficiently smooth velocity 
field v(r). In a smooth velocity field, an initial singularity 
of the generalized vorticity could not disappear in subsequent 
moment of time, it could be only transported and deformated 
by the fiuid flow. As far as the equations of motion for an in- 
viscid substance are time-reversible, it follows from this, that 
the formation of singularity from a smooth initial data also 
would not possible. So, the existence and possible types of 
collapse of vortex lines, in the systems with the quadratic on 
ft Hamiltonians, depend on the asymptotics of the interaction 
function G(ri,r2) at ri ^ r2. So, for a better understanding 
of the collapse problem in hydrodynamics, it is reasonable to 
investigate such systems, for which the interaction function 
has the asymptotics G ~ |ri — r2|^'' , and the exponent q is 
not necessary equal to the unity. 

What is the influence of viscosity on the structure of col- 
lapsing region? How does the given type of collapse have 
effect on turbulent spectra? That are only few of the most 
important questions, which need investigation. It would be 
interesting to verify numerically our hypothesis about possi- 
bility of (quasi-)inertial collapse in ideal hydrodynamics, both 
in Eulerian variables and in the vortex line representation. 
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